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The fluctuation in electric current in nonequilibrium steady states is investigated by molecular 
dynamics simulation of macroscopically uniform conductors. At low frequencies, appropriate 
decomposition of the spectral intensity of current into thermal and excess fluctuations provides a 
simple picture of excess fluctuations behaving as shot noise. This indicates that the fluctuation- 
dissipation relation may be violated in a universal manner by the appearance of shot noise for 
a wide range of systems with particle or momentum transport. 
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In equilibrium states, the fluctuation of an observable 
is related universally to a linear response function by the 
fluctuation-dissipation relation (FDR).^ In nonequilib- 
rium steady states (NESSs), the FDR is often violated, 
and 'excess fluctuation' (XF) appears. XF plays cru- 
cial roles in many flelds of physics, including single elec- 
tron tunneling,^ the squeezing of photons,'^' ^ the mea- 
surement of fractional charge,^ and the determination of 
the fundamental limits of quantum interference devices.^ 
However, unlike equilibrium fluctuation, it is not yet well 
understood whether universal properties exist in XF.^^ 

Experimentally, the FDR violation is hardly observ- 
able in heat conduction because a convection current or 
a phase transition is induced for large temperature differ- 
ence (which drives heat conduction) before XF becomes 
detectable. In contrast, the violation has been widely ob- 
served in systems with particle (or momentum) trans- 
port, such as electric conductors and photoemitting de- 
vices. '^"^^ We therefore consider such systems. 

Among such systems are simple systems, includ- 
ing mesoscopic conductors, ^"^^ conductors with junc- 
tions^^' ^"^ (e.g., tunnel and PN junctions), and light- 
emitting diodes.^' ^ These systems are simple in the sense 
that the number of electron modes is small and/or many- 
body interactions are unimportant and/or dissipation is 
negligible and/or the principal origin of XF is localized in 
certain mesoscopic regions. XF generated in such a case 
takes the form of shot noise. ^"^^ Here, the term 'shot 
noise' is used in a wide sense, which stands for fluctu- 
ation whose spectral intensity Si is proportional to the 
absolute value of average flux, |(/)|.^^ The ratio W of Si 
to its Poissonian value is called the Fano factor, which 
takes various values depending on the details of the sys- 
tems. ^"^"^ 

The situation is completely different for uniform 
macroscopic conductors, for which the assumptions made 
in refs. 3-17 do not hold. Although the FDR violation is 
hardly observable in uniform metals, it is widely observed 
in uniform semiconductors.^^ Most experiments on the 
latter showed that XF is dominated by 1// noise, which 
is proportional to {I)"^ }^ Although shot noise may also 
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exist in such systems, it would be masked by 1// noise^*' 
because the latter increases more rapidly with increasing 
|(/)|. However, the origin of 1// noise is believed to be 
imperfections in samples, such as the ffuctuation in car- 
rier number and the migration of impurities, which result 
in a strong sample dependence of 1// noise. Since im- 
perfections in samples are of secondary interest in funda- 
mental physics (nonequilibrium statistical mechanics), a 
natural question arises: What fluctuation appears in per- 
fect samples? In this paper, we address this question and 
report a property of XF that may be universal. 

The models and results of the previous works on meso- 
scopic conductors^"^^ are not applicable to macroscopic 
conductors, because, as mentioned above, many assump- 
tions that do not hold in macroscopic conductors have 
been made in those works. We therefore take a different 
approach. That is, we use molecular dynamics (MD) sim- 
ulation on a model that we believe captures the essential 
elements of macroscopic conductors. This enables us 
to study the NESSs of perfect samples, without making 
the assumptions made in the works on mesoscopic con- 
ductors. Since we can vary the values of the parameters 
to a great extent, we are able to present results that may 
indicate a universal character. 

Except at low temperatures, quantum effects seem to 
play minor roles in macroscopic conductors far from equi- 
hbrium, because of the strong decoherence. We there- 
fore use the classical model of electric conduction pro- 
posed in ref. 21, which describes doped semiconductors 
at room temperature well.^^ The system includes three 
types of classical particles, which we call electrons (each 
with mass rUe and charge e), phonons (each with mass 
TTip), and impurities. Their number densities are denoted 
by ric. Up, and rii, respectively. For simplicity, we assume 
a two-dimensional system, the size of which is Lx x Ly. 
In the x-direction, we apply an external electric field act- 
ing only on electrons, and impose the periodic boundary 
condition. The boundaries in the y-direction are poten- 
tial walls for electrons and thermal walls for phonons. 
The thermal walls reflect phonons with random velocities 
sampled from an equilibrium distribution with tempera- 
ture Tq?^ This enables phonons to carry heat constantly 
out of the system thereby keeping the system in a NESS. 



2 J. Phys. Soc. Jpn. 



Letter 



Author Name 



Impurities are immobile and play a role of providing ran- 
dom potential. 

We assume short-range interactions among all parti- 
cles. Since interaction potential is well characterized by 
its scattering cross section, detailed forms of the poten- 
tial are expected to be irrelevant when studying general 
nonequilibrium properties. Therefore, we here take a sim- 
ple form, fco(max{0, cij/})^/^. Here, ko is a constant and 
dji = Rj + Ri — \rj — ri\ is the overlap of the poten- 
tial ranges. Rj is the radius of the potential range {R^, 
Rp, and Ri for an electron, phonon, and impurity, re- 
spectively), and rj is the position of the j-th particle. 
We can change the strength of scattering by varying Rj 
(and particle number density). 

This model corresponds to a perfect sample because 
the total number of carriers does not change and because 
impurities do not move. This system is macroscopically 
uniform although the translational invariance is broken 
by impurities and the thermal walls for phonons. Fur- 
thermore, the model and results are also applicable to 
systems that have a mass flow of neutral particles. 

Wc use units in which irie, Re, e, the Boltzmann con- 
stant, and a reference energy are unity. Regarding the 
other parameters, the main result, eq. (5), is insensitive 
to their values, as will be shown later. We here fix Rp = 1, 
To = 1, and ko = 4000; the other parameters are varied 
to illustrate the possible universality of the result. 

To the investigate nonequilibrium states of this model, 
we perform MD simulation using Gear's fifth-order 
predictor-corrector method. The time-step width is set 
to 10~^. The initial position of each particle is ran- 
domly arranged so as not to be in contact with the other 
particles, and the initial velocities of the electrons and 
phonons are given by the Maxwell distribution with tem- 
perature Tq. We calculate various quantities after the 
system reaches a NESS. 

The electric field applied to the system is composed 
of a time-independent field E, which is varied in a 
wide range, and a time-dependent field ef{t), which is 
small. The electric field induces electric current I{t) = 
encLyV^{t), where is the velocity in the x-direction 
(i.e., along the electric field) of the center of mass of 
electrons. We take e 7^ only when we calculate the dif- 
ferential response function /x(t — t; E) of a NESS, which 
is defined by 

{SI{t))E,e= I dr ii{t-T;E)L,ef{T) + 0{e^), (1) 
J —00 

for f>T and by fj,{t — t;E) = for f < r. Here, 51 = 
I — (/) E,o, and (• • • ) E,e denotes the average at the NESS 
in the electric field E + sf{t). The convolution theorem 
yields fl{uj; E) = \v[a^^t){5I{uj))E,sl LxSf^oj), where the 
tilde denotes the Fourier transform. Note that jl{u);E) 
differs from that in an equilibrium state, fi{u);0). 

Wc are mainly interested in the current fluctuation 
that is characterized by the spectral intensity Si{lu;E) 
of I{t) for e = 0. By the Wiener-Khinchine theorem,^ 
Si{uj;E) is equal to the Fourier transform of the auto- 
correlation function {6I{t)SI{0)) e,o of current. In equilib- 
rium states {E = 0), the FDR, Si{lj;0) = 2TRe/i(a;; 0), 
holds for all w.^ Here, T is the temperature of the con- 
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Fig. 1. (a) Si{uj;E] and 2ToRcft(u); E) for E = 0. The inset 
shows {I)e,o versus E. (b) Si(uj;E), 2Te{E)Rejl{w; E), and 
2ToKejj,{(jj; E) for E = 0.06 (nonhnear response regime). In 
these simulations, rup = 1, Ri = 0.5, Lx = 750, Ly = 125, 
Ue = Up = 0.016, and ni = 2/375. The data points are the aver- 
ages of five samples (impurity configurations) and the error bars 
are the standard deviations among them. 

ductor, which is equal to To when E = 0. We plot both 
sides of this relation in Fig. 1(a), and confirm that it 
holds in our simulation. 

When larger E (7^ 0) is applied, (/) e.q becomes non- 
linear with E, as shown in the inset of Fig. 1(a). In such 
NESSs, we find that the FDR is violated, i.e., for any T 
that is independent of co, 

Si{u); E) 7^ 2TRe/},(w; E) for some w. (2) 

This is demonstrated in Fig. 1(b), which shows 6*7(0'; E), 
2ToReil{oj; E), and 2Te{E)Reil{oj; E) in a nonlinear re- 
sponse regime. Here, Te{E) = me(« - {Ve)E,of)E,o 
is a kinetic temperature of electrons (v^ is the veloc- 
ity of an electron in the cc-direction). When we employ 
2ToRe/i(w; E) as the right-hand side (RHS) of the FDR, 
the violation of the FDR is observed in a wide frequency 
range. When we use 2Te{E)Rcjl{uj; E) as the RHS, the 
violation is observed at low frequencies (w <IC wo)^^ while 
the RHS coincides with Si{u!;E) at higher frequencies 
(uj ^ luq), where is the crossover frequency between 
the regimes of FDR violation and validation. These data 
also show that the FDR is violated for any definitions of 
T that is independent of lo. 

Now we discuss the main finding of this paper. Since 
we have seen that the FDR violation is manifested at 
lower frequencies, wc look at the low-frequency region 
(w <^ uJo)- Among many possible definitions of 'thermal 
fiuctuation' of / for i? 7^ 0, we employ 

StiXLu;E)=2ToReil{Lo;E), (3) 

which is the RHS of eq. (2) with T = Tq. Using this, we 
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Fig. 2. Excess fluctuation 5exs at a low frequency, plotted against 
{I)e,o- The dotted line represents the asymptote, W(|(/)b,o| — 
7o), fitted with the four data points at larger values of {I)e,o- 
The parameters of the simulation and the meaning of the error 
bars are the same as those in Fig. 1. 
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decompose the total fluctuation Sj into two parts: 

5/(0;; E) = Sth{oj; E) + 5exs(w; E). (4) 

Since the thus-defined 5exs quantifies the FDR violation, 
wc call it excess fluctuation. In Fig. 2, we plot 5*0X3 for 
LO ~ 0.002 as a function of {I)e,o- [We can translate a 
function of E into a function of (I) e,o because of the one- 
to-one correspondence between E and {I)E.a-] Since the 
FDR holds in equilibrium states, ^exs — when {I)e,o 
is small. As {I)e,o increases, Sq^s exhibits a crossover 
behavior from near equilibrium to far from equilibrium 
as 

^exs-| w{\{I)e,o\-Io) (K/)b,o|»/o), ^' 

where Iq is a certain crossover vahie of the current. In 
the latter region {\{I)e,o\ ^ -^o), •S'exs takes the form of 
shot noise, where W is the Fano factor. 

We have thus found that the dominant mechanism that 
breaks the FDR is the appearance of shot noise. To con- 
firm that this observation holds widely for the model con- 
sidered here, we also study Soxs in the following cases: (i) 
another impurity density, ni = 0.016, (ii) other linear di- 
mensions Lx (along E) of the system, = 375, 300, 
187.5, and 150, and (iii) the values of the other parame- 
ters are changed significantly (e.g., Ue = 0.008, rUp = 10, 
and iii = 2). (iv) The thermal walls for phonons are set 
away from the boundaries for electrons, as shown in the 
top- left inset of Fig. 3(b). 

Figure 3 shows the results in case (iv). In this case, 
the local phonon temperature Tp around the boundaries 
for electrons is markedly different from Tq, as shown in 
Fig. 3(a), where Tp = mp{{v^-{v^)Efif)Efi {Vp is the lo- 
cal phonon velocity in the x-direction). Despite this fact, 
5*6X8 is well-fitted again by eq. (5), as shown in Fig. 3(b), 
if we define thermal fiuctuation again by eq. (3) using Tq. 
Furthermore, we have found, although the data are not 



Fig. 3. (a) Local phonon temperature Tp and (b) excess fluctu- 
ation Sexs for LU ~ 0.002, plotted against {I)e,o, for a system 
where the thermal walls (at y = and 300) for phonons are 
set away from the potential walls (at y = 87.5 and 212.5, dash- 
dotted lines) for electrons, as shown in the top-left inset of (b). 
In (a), the solid lines from bottom to top correspond to the data 
for E=0, 0.01, 0.02, 0.06 and 0.12. In (b), the dotted fine repre- 
sents the asymptote W e.q\ ~ ^o)- The meaning of the error 
bars is the same as that of the error bars in Fig. 1. Bottom-right 
inset: {I)e,o versus E for this system. We take nip = 1, i?i = 0.5, 
= 375, tie = 0.016, Up = 19/1125, and ni = 2/375. 

shown here, that eq. (5) also holds well in cases (i)-(iii) 
(except when the densities are so high that a liquid-solid 
phase transition takes place). Note in particular that the 
validity of eq. (5) in case (iii) suggests that it holds in- 
dependently of details of the models, because case (iii) 
naturally includes, for example, the case where the Rj 
are specific functions of Ue- 

The above observations strongly indicate the robust- 
ness of eq. (5). Note that this possible universality is 
visible only when thermal fluctuation in nonequilibrium 
states is appropriately defined as eq. (3). In fact, we have 
found (although the data are not shown here) that the 
possible universality is obscured if we use Te{E) instead 
of To in thermal fluctuation. 

Using the results in case (ii), we also investigate the 
Lx dependences of W and Iq. We evaluate W and /q 
by fltting the numerical results of 5cxs for large |(/)b,o| 
to the asymptotic form of eq. (5). In Fig. 4, we show 
WLx and Iq versus L^- We see that WL^ is almost in- 
dependent of Lx, i.e., W ~ 1/Lj;. This agrees with the 
partial result for macroscopic conductors in ref. 9 (how- 
ever, see note^^), and coincides with the results for long 
mesoscopic conductors.'''* Furthermore, we observe that 
Iq is almost independent of L^, although the error bars 
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Fig. 4. WLx (circles; left axis) and Iq (squares; right axis) for 
several sizes of the system. The particle densities are the same 
as those in Fig. 1. 



are somewhat large. 

By conibining the present results with the results on 
simple systems, ^"^'^ we conjecture that the FDR is vio- 
lated not in a random and system-dependent manner but 
in a universal manner by the appearance of shot noise, for 
a wide range of systems from mesoscopic to macroscopic. 
All details of individual systems are absorbed into W, 
Iq, and the differential response function Rejl (by which 
thermal fluctuation is deflned). The origin of current fluc- 
tuation in the present model is the chaotic behavior of in- 
teracting many particles in classical systems, while that 
in the simple systems^"^^ is essentially the probabilis- 
tic nature of quantum or thermal-activation processes of 
noninteracting particles. Despite such a marked differ- 
ence, Sexs takes an identical form in all systems. This 
observation may be used as a touchstone in nonequi- 
librimn thermodynamics or statistical mechanics beyond 
the linear response theory. 

Note that the present results could never be obtained 
by a naive perturbation expansion, in powers of the driv- 
ing force E, about an equilibrium state. For example, 
the relation 5exs oc \{I)\ suggests that such a power 
series would not converge for large E. Using MD sim- 
ulation, we have successfully investigated such a 'non- 
perturbative regime.' Our results may be confirmed ex- 
perimentally, for example, in high-quality doped semi- 
conductors, which may be prepared by modulation dop- 
ing, at room temperature. 

In conclusion, we have presented a study of excess fluc- 
tuations in a nonequilibrium system and found that the 
fluctuation-dissipation relation is violated in a manner 
that may be universal. We hope that our work will stim- 
ulate further research that will test the correctness of 
this conjecture for wider classes of systems. 
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